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Abstract 

We present a general framework for analyzing spatially inhomogeneous cosmological dynamics. 
It employs Hubble-normalized scale-invariant variables which are defined within the orthonormal 
frame formalism, and leads to the formulation of Einstein's field equations with a perfect fluid matter 
source as an autonomous system of evolution equations and constraints. This framework incorporates 
spatially homogeneous dynamics in a natural way as a special case, thereby placing earlier work on 
spatially homogeneous cosmology in a broader context, and allows us to draw on experience gained 
in that field using dynamical systems methods. One of our goals is to provide a precise formulation 
of the approach to the spacelike initial singularity in cosmological models, described heuristically by 
Belinskii, Khalatnikov and Lifshitz. Specifically, we construct an invariant set which we conjecture 
forms the local past attractor for the evolution equations. We anticipate that this new formulation 
will provide the basis for proving rigorous theorems concerning the asymptotic behavior of spatially 
inhomogeneous cosmological models. 
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1 Introduction 



Scales and scale invariance play a crucial role in practically all branches of physics, and general relativity 
(GR) and cosmology are no exceptions.^ In these cases one is interested in self-gravitating systems, 
which in the cosmological context requires a matter model as well as a spacetime description. This in 
turn requires consideration of scales. In modern cosmology one assumes that (today) there exists a global 
scale — that of the particle horizon. The empirical data are usually interpreted as follows: on sufhciently 
large spatial scales, say a few percent of the particle horizon, everything looks statistically roughly the 
same in all directions. Combined with the Copernican principle ("we are not located at a special place"), 
this suggests that one can replace a very complicated matter distribution by a simple one: a smooth 
distribution that is spatially homogeneous and isotropic, obtained by averaging over sufficiently large 
spatial scales. Then it is further assumed that one can also approximate the geometry of the spacetime 
by a spatially homogeneous and isotropic geometry, i.e., one assumes that the geometrical features trace 
those of the matter and that possible "excited geometrical modes" , like gravitational waves, are negligible 
on these scales. This then leads to modeling the cosmological spacetime by a Robertson- Walker (RW) 
geometry. 

'Electronic address: Claes .UgglaSkau. se 
t Electronic address: H. van. Elst{Bqmul . ac .uk 
-^Electronic address: jwainwriOmath.uwaterloo . ca 
^Electronic address: ellisamaths .uct . ac .za 

^See, e.g., the recent Resource Letter by Wiesenfeld on scale invariance in physics and beyond |37| . 
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The assumption of a RW geometry subsequently forces the summed matter content to take the form 

of a perfect fluid through Einstein's field equations (EFE) . Although the mathematically simplest matter 
model is a single perfect fluid, we need more complex matter models to describe the real Universe. Indeed, 
matter in the Universe consists of many components: at least (i) radiation (photons), (ii) baryonic matter, 
(iii) neutrinos, (iv) dark matter, and (v) dark energy/quintessence (other components like cosmological 
magnetic fields are usually neglected). Once the matter content has been specified and equations of state, 
scalar field potentials, particle distribution functions, etc., have been chosen, the evolution of the model 
is determined by the EFE and the matter equations, e.g., the evolution equation for a scalar field. This 
then leads to a Friedmann-Lemaitre (FL) model for the Universe. 

The next step, aimed at describing the actual inhomogeneous Universe, is to perturb the FL model 
and describe the evolution of large-scale structures in the Universe, which appear at many scales — 
filaments and voids, superclusters of galaxies, galaxies etc. But it is generally believed that linear per- 
turbation theory can account for them all on these large scales. And, indeed, the FL scenario and the 
linear perturbations thereof ( "almost-FL models" ) form a remarkably successful framework — it seems to 
consistently account for present observational evidence, at least over sufficiently large smoothing scales. 
Moreover, it forms an interpretational framework that encourages and steers further observations. These 
are currently focused on determining the various density contributions Qi (including J^a for the cosmo- 
logical constant), and the spectrum and growth of density perturbations. This is the simplest scenario 
consistent with current observations. 

Nevertheless, there are issues that need elucidation that by necessity lie outside the domain of the 
standard almost-FL picture. Here are some of them; 

• To investigate the constraints observations impose on the spacetime geometry of the Universe 
requires investigating a hierarchy of more general models, perhaps characterized by assumed "priors" 
(where removing a prior necessarily involves looking beyond ones favorite model, with the hope of 
getting further support for it!). 

• To understand how special the FL models are and put them into a broader context requires looking 

beyond them. 

• While the Universe is well described by almost-FL models at present, this may not always have 
been true, and it may not remain true in the far future; in particular, we would like to know the 
largest class of models that can look like a FL model at some stage of their history. 

• Using a FL and a linearized FL scenario rules out possible non- linear effects, but these may domi- 
nate; e.g., as structures become much more dense while aggregating to smaller scales. 

• GR is a highly non-linear theory; to prove that the linear theory is correct requires going beyond 
linear perturbations. 

• The averaging and fitting procedures motivating the FL models do not a priori commute with the 
EFE, i.e., starting with an inhomogeneous model and smoothing it does not necessarily lead to the 

model that has been smoothed from the outset and then perturbed. This gives rise to a number of 
questions, e.g., can inhomogcneitics affect the overall evolution? How do they affect observations? 

• There are deep connections between GR, cosmology and thermodynamics, e.g., relating (gravita- 
tional) entropy and the arrow of time. To better understand such connections requires a state space 
picture describing the set of solutions, where one can examine coarse-graining and existence of at- 
tractors on this state space, toward which the evolving cosmological models move. Since entropy 
requires counting of possible states this requires looking at models beyond FL. 

• What is the detailed nature of possible singularities? A better understanding of generic features 
of singularities and their dependence on matter content and initial data might shed light on how 
the real Universe evolved initially. There might also exist at least a local mathematical connection 
between the initial singularity and the singularities of gravitational collapse. To understand such a 
relationship, or its non-existence, again requires an inspection of cosmological solutions beyond the 
restrictions imposed by RW geometries. 
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• A better classical understanding of singularities might help to produce gravitational theories with 
greater domains of validity; e.g., finding asymptotic symmetries of the field equations when ap- 
proaching singularities may provide sufficient structure to asymptotically quantize the theory in a 
regime where quantum gravity is supposed to be of importance. 

Thus, there is ample motivation to probe a larger subset of the cosmological solution space of the EFE 
than just the almost-FL models. Our first goal in this paper is to develop a framework for this purpose. 

In view of the above-mentioned importance of scale invariance in physics, we propose to introduce 
scale-invariant variables, and to describe the evolution of a cosmological model by an orbit in an infinite- 
dimensional dynamical state space, governed by first-order (in time) autonomous evolution equations 
derived from the EFE and the matter equations. The behavior of the model in the asymptotic regimes, 
i.e., near the initial singularity and at late times, can then possibly be described by a past attractor and 
a future attractor of the evolution equations. 

The appropriate mathematical vehicle for implementing this proposal is the orthonormal frame for- 
malism, since (i) it describes the essential degrees of freedom of the gravitational field in a coordinate 
independent manner, (ii) orthonormal frame vectors provide local reference scales and so it allows one to 
naturally introduce scale-invariant variables, and (iii) it leads directly to first-order (in time) autonomous 
evolution equations. 

To be more specific, by a cosmological model we mean a four-dimensional spacetime manifold M 
endowed with a Lorcntzian metric g which satisfies the EFE with an appropriate matter/energy distri- 
bution. We assume the existence of a local foliation of the spacetime manifold by a one-parameter family 
of spacelike 3-surfaces with a future-directed unit normal congruence u. We naturally choose this unit 
vector field to be the timelike vector field in the orthonormal frame. We assume that the cosmological 
model is expanding, i.e., the volume expansion rate'^ Q of the normal congruence is positive. Because we 
are working in a cosmological setting we will replace Q by the Hubble scalar^ H — ^Q. 

In the orthonormal frame formalism, the frame vector components and the commutation functions are 
the basic dynamical variables for describing gravitational fields, and they each have physical dimension* 
[ length ]^^. The Hubble scalar also has physical dimension [length]"^, and constitutes the natural 
cosmological length scale through the Hubble radius H^^. This fact motivates one of the key steps 
in our approach, namely, the introduction of Hubble-normalized variables by dividing the frame vector 
components and the commutation functions by H. Curvature quantities such as the matter density and 
the orthonormal frame components of the Weyl curvature tensor have physical dimension [ length ]^^, 
and hence are normalized by dividing by H^. This process of Hubble-normalization has two important 
consequences. Firstly, dimensional variables are replaced by dimensionless ones, leaving H as the only 
variable carrying physical dimensions. Secondly, one is essentially factoring out the overall expansion 
of the Universe, thereby measuring the dynamical importance of physical quantities (e.g., the matter 
density) relative to the overall expansion (cf. Kristian and Sachs HHI)- This choice also provides a link 
between mathematical analysis and observation, since key observational variables are Hubble-normalized. 
Earlier investigations of the asymptotic dynamics of cosmological models using scale-invariant variables 
dealt with spatially homogeneous (SH) cosmologies, i.e., models that admit a three-parameter group of 
isometrics acting on spacelike 3-surfaces (see Wainwright and Ellis (WE) ^-nd other references therein) , 
and the so-called G2 cosmologies, which admit a two-parameter group of isometrics acting on spacelike 
2-surfaces (see Refs. ^Kj and |23j)- The framework that we develop in this paper generalizes a program 
that has been extremely successful in a SH context to a completely general spatially inhomogeneous 
setting, i.e., to models that admit no isometries. Indeed, the SH and G2 cosmologies will be incorporated 
in a natural way as invariant sets of the infinite-dimensional Hubble-normalized state space. For brevity, 
and to emphasize that they admit no isometries, we shall refer to the models under consideration as Go 
cosmologies. 

The framework that we are developing can be used to study both asymptotic regimes in an ever- 
expanding cosmological model. In this paper we focus on the initial singularity. The problem of asymp- 
totics in spacetimcs that exhibit no symmetries poses a formidable challenge. Nevertheless, concerning 

■^Here = — (trfc), where © is the volume expansion rate of the normal congruence and (trfc) is the trace of the extrinsic 
curvature of the spacelike 3-surfaces. 

^When evaluated at the present epoch, the Hubble scalar equals the Hubble constant Hq, familiar from observational 
cosmology. 

*We will use units such that Newton's gravitational constant G and the speed of light in vacuum c are given by H-kG/c^ = 1 
and c = 1. 
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the existence of singularities, remarkable progress was made several decades ago by Penrose and Hawk- 
ing, leading to their singularity theorems |21[I22| . However, the singularity theorems do not tell us much 
about the nature of the singularities. Detailed asymptotic analysis, using the full EFE, is required for 
this purpose. To date rigorous results have, with few exceptions, been confined to cosmological mod- 
els with isometrics, in particular SH and G2 cosmologies. We shall discuss these results in Sec. As 
regards initial singularities in Go cosmologies, heuristic results were obtained by Belinskii, Khalatnikov 
and Lifshitz (BKL) [32] by making ad hoc metric assumptions that were subsequently inserted into the 
EFE with the purpose of showing that they were consistent. This analysis led to a remarkable, although 
heuristic, conjecture that has become part of the folklore of relativistic cosmology. 

The BKL conjecture: 

For almost all cosmological solutions of Einstein's field equations, a spacelike initial singularity is 
vacuum- dominated, local and oscillatory . 

For cosmological models with a perfect fluid matter source, the phrase "vacuum-dominated" , or, 
equivalently, "matter is not dynamically significant" , is taken to mean that the Hubble-normalized matter 
density (i.e., the density parameter f2) tends to zero at the initial singularity. The phrase "for almost all" 
is needed because there are a number of exceptional cases. Firstly, if the perfect fluid has a stiff equation 
of state, the density parameter does not tend to zero (see Andersson and Rendall Secondly, there 
is a special type of initial singularity called an isotropic initial singularity, in the neighborhood of which 
the solution is approximated locally by a spatially flat FL model (see Goode and Wainwright Q^), with 
the result that the density parameter tends to the value 1. Isotropic initial singularities, however, only 
arise from initial data that form a set of measure zero.^ 

The word "local" in the BKL conjecture means heuristically that the evolution at different spatial 
points effectively decouples as the initial singularity is approached, with the result that geometrical 
information propagation is asymptotically eliminated. It is natural to describe this phenomenon as 
asymptotic silence of the gravitational field dynamics. We shall refer to the associated initial singularity 
as being a silent initial singularity.^ 

The word "oscillatory" in the BKL conjecture means that the evolution into the past along a typical 
timeline passes through an infinite sequence of Kasner states, generalizing the behavior first encountered 
in the so-called Mixmaster models (SH models of Bianchi Type-IX; see Misner 

Our second main goal in this paper is to give a precise statement of the BKL conjecture, within the 
framework of the Hubble-normalized state space. 

The plan of this paper is as follows. In Sec. [21 we derive the Hubble-normalized evolution equations 
and constraints for Gq cosmologies that arise from the EFE and the matter equations. In Sec.[2|we make 
a choice of gauge and then describe some features of the Hubble-normalized state space, in particular 
the SH invariant set and the silent boundary. We then define the notion of a silent initial singularity. 
In Sec. 21 by analyzing the dynamics on the silent boundary, we are led to construct an invariant set 
which we conjecture is the local past attractor for Go cosmologies with a silent initial singularity, thereby 
making precise the notion of an oscillatory singularity. In Sec.[5lwe consider various classes of cosmological 
models with isometrics and use the past attractor to predict whether the initial singularity is oscillatory 
or not. We conclude in Sec. [^1 with a discussion of silent initial singularities and the BKL conjecture, and 
raise some issues for future study. Useful mathematical relations, such as the propagation laws for the 
constraints and expressions for the Hubble-normalized components of the Weyl curvature tensor, have 
been gathered in an appendix. 

2 Evolution equations and constraints 

In this paper, we consider spatially inhomogeneous cosmological models with a positive cosmological 
constant. A, and a perfect fluid matter source with a linear barotropic equation of state. We thus have 

p(/i) = (7-l)A, (1) 

^ There is in fact a wide variety of known special SH and G2 solutions in which the initial singularity is matter-dominated, 
i.e., n does not tend to zero. Like the isotropic initial singularities, these singularities only arise from initial data that form 
a set of measure zero. 

^For further discussions on the suppression of information propagation and asymptotic silence, see Ref. 1151 . 
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where fl is the total energy density (assumed to be non-negative) and p the isotropic pressure, in the rest 
3-spaces associated with the fluid 4-velocity vector field m, while 7 is a constant parameter. The range 

1 < 7 < 2 

is of particular physical interest, since it ensures that the perfect fluid satisfies the dominant and strong 
energy conditions and the causality requirement that the speed of sound should be less than or equal 
to that of light. The values 7=1 and 7 = | correspond to incoherent pressure-free matter ("dust") 
and incoherent radiation, respectively. In cosmology it is natural to single out a future-directed timelike 
reference congruence eo = m of unit magnitude. This gives rise to a (H-3)-decomposition of the perfect 
fluid energy-momentum-stress tensor 

Tab = Ufc + 2 g(a Ufo) + P /lab + TTab , (2) 

with 

M = r^G+A P = g;M(7-i) + (i-|7)^'']m 

(3) 

= jG^^fiv" TTab = 7G_^V^^(aWb> 

The vector field v, which represents the peculiar velocity of the fluid relative to the rest 3-spaces of Bq, 
is defined by 

M° :=r(u"+w") , VaU^^O, (4) 

with the Lorentz factor given by 

r := -j=L= , vav'' . (5) 

V 1 — f 

The scalars G± (we shall require G_ later) are defined by 

G± :=l±(7-l)i'^ (6) 

To obtain an orthonormal frame, { }a=o,i.2,3j we supplement the timelike reference congruence Bq with 
an orthonormal spatial frame { }q=i.2,3 in the rest 3-spaces of eg. The frame metric is then given 
by '7ab = diag [ — 1, 1, 1, 1 ]. In the orthonormal frame formalism, introduced in relativistic cosmology, 
among others, by Ellis the basic variables are the frame vector components, the commutation 

functions associated with the frame and the matter variables, and the dynamical equations are provided 
by the EFE, the Jacobi identities and the contracted Bianchi identities (the latter, for a perfect fluid, 
corresponding to the relativistic extension of Euler's equations). We will make use of an extended version 
of this formalism given by van Elst and Uggla ^^l- The dynamical equations consist of two sets, those 
containing the temporal frame derivative eo, which we refer to as evolution equations, and those not 
containing Cq, which we refer to as constraints. 

To convert the dynamical equations of the orthonormal frame formalism to a system of partial differ- 
ential equations (PDE), it is necessary to introduce a set of local coordinates {a;^}^-o,i,2,3 = {t, a;*}i=i,2,3- 
We do so by adopting the standard (3-|-l)-approach (see, e.g., Refs. 38, and ^Sf^). Here Bq is assumed to 
be vorticity-free and, thus, hypersurface-orthogonal. As is well-known, this gives rise to a local foliation 
of the spacetime manifold by a one-parameter family of spacelike 3-surfaces, S:{t = constant}. The 
(3-f l)-approach leads to the following coordinate expressions for the frame vector fields (cf. Ref. 53)- 

Bo = N-^ {dt - di) , e„ = , (7) 

where N and iV' are known as the lapse function and the shift vector field, respectively. 

2.1 Dimensional equation system 

We now present the dynamical equations as given in Ref. simplified by the assumption that Bq is 
vorticity-free {uj" — 0). We begin with the commutator equations, which serve to introduce the basic 
gravitational field variables, and which will later be used to derive some additional evolution equations 
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and constraints. 
Commutator equations: 

[eo,e„](/) = Ua.eo{f)~{HSj + aj -e^/n'')epif) (8) 
= (Ccom)a/3(/) := [e„,e^](/)-(2a[„,50f + e„;55n^T)e^(/) , (9) 

with / denoting an arbitrary real- valued spacetime scalar. Here H is the Hubble scalar which is related to 
the volume expansion rate 8 of Cq according to H := ^Q. The quantities and are the acceleration 
and shear rate of Cq, respectively, while fl" describes the angular velocity of the spatial frame { } along 
the integral curves of Bq relative to a Fermi-propagated one. The quantities a" and determine the 
connection on the spacelike 3-surfaces S:{t = constant}. 

Einstein's field equations, Jacobi identities and contracted Bianchi identities (Euler's equations): 
Evolution equations: 

eo{H) = i(^^^f^"/3)_ i(^ + 3p)+ 1 A+i(e„+7i„-2a„)(u") (10) 

eo(a") = -(i/(5"^i + a"/3-e%/3f^^)a''-i(e^ + W/3)(2i?(5"'3-(T"'5-e"'5^r!^) (11) 

+ e'"'^" [ (e^ + u^~ 2a^) {n'^\) + 2fl^ a'^\]+ tt"^ + (S''^" + + a<") (ti''^) (12) 
eo(n"'^) = -{H6^"s-2<7^°'s-2e^s^°'^',)n'^''^ -{e^+u^)ie'-'^^°'a^^'>s-S''^°'i}'^'> +6'''^n^) (13) 
eo{^i) - -3i/(/i + p)-(e„ + 27i„-2a„)(g")-(a„/3^"'') (14) 

eoiv'') = [-7«"eo(M) + (G-5"^ + 2(7-l)t;"«0)eo(g'')] , (15) 

7 G_ fj. 



where 



- {ep + up~ iap) (^"^5) + e'^f"^ nps tt/ . (16) 



Constraints: 



= (Cg) := 2(2e„-3a„)(a")-(n„;3'^"^) + i(n„")2 + 6//2_(^^^^"/3)_2M-2A (17) 

= {CcT — -e^(2i7^"'3_^a/3)_3^^^a/3_ga/37„^^^^5^ga (^g) 

= (Cj)" e;3(«"'' + e"^'^a7)-2a^n"'3 . (19) 

Note that we are not provided with evolution equations for any of the'^ coordinate gauge source functions 
N and iV' (which reside in eg) or the frame gauge source functions and fi". Note also that these ten 
gauge source functions do not appear in the constraints. Independent of a choice of gauge (to be discussed 
in Sec. inj, the evolution equations ((HJ and HlOfl (|16|l propagate the constraints ^ and (|17(l - (|19|l along 
the integral curves of Cq according to Eqs. (|139|I - H142|I in the appendix. 

There are, in addition, two gauge constraints that restrict four of the gauge source functions, given 

by 

= {C^r ■■= [e"''Me/3-a/3)-n"^]u^ (20) 
= (aOa N-^e^{N)-u^ . (21) 

The former is a consequence of assuming eg to be vorticity-frec, the latter follows from Eq. Q upon 
substitution of Eq. ((T)). The propagation of the gauge constraints along the integral curves of eg can 
be established once a choice of temporal gauge has been made (as this determines what the currently 
unknown frame derivatives eo(u") and eo{N) should be). 

'^Employing the terminology of Friedrich 1161 . Sec. 5.2. 
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2.2 Scale-invariant equation system 

We now introduce Hubble-normalized frame, connection and matter variables as follows: 

^0 := eo , da ■= J^^a , (22) 

{U",Eap, A'', Na0. R"} := {u" , a^fi, a" , n^p, n'^ }/ H (23) 
{ n, n^, P, Q", Ua0} := { A, p, g", 7rap}/{3H'~) . (24) 



It follows from Eq. © that 



P = G;M (7-l) + (l-|7)«' , Q" = -^f7z;", Uap^^nvi^avp) ■ (25) 

Expressing the Hubble-normalized frame derivatives do and da with respect to the local coordinates 
introduced in Eq. Q leads to: 

do = AA-i {dt - d,) , da = Ej , (26) 

where 

Af:=NH, Ea' :- ^ . (27) 

rl 

In order to write the dimensional equation system in Hubble-normalized form, it is necessary to introduce 
the deceleration parameter q and the spatial Hubble gradient Tq, defined by 

(q + l) := ~^doH (28) 

ra -^daH. (29) 

The definition H28|) . together with Raychaudhuri's equation (|10|l and Eqs. H23|l and (|24|l . lead to the 
following key expression for q: 

q = 2Y? + \{^ + ?,P)-nA- \ [da -ra + ila- 2A„) U" , (30) 
where := i (S^.^E"'')^ 

We now use Eqs. (|28|l and 129|1 to write the commutator equations © and in Hubble-normalized 
form. The result is 

[do,da]f = -{ra^Ua)dof+{q5a^~^a^ + ea^^R^)dpf (31) 
= (Ccondcpif) ■■= [da,dp]f-[2{r[a+A[a)Sp]''+eai3sN''^]d^f . (32) 

We now write the evolution equations (|10|) - (|16|l and the constraints p7|) - (|19|l in Hubble-normalized 
form.® 

Evolution equations: 

aoA" - (gr^-S^^-f e%0i?^)A'5-i(a/3-r^-f [/^)(2,5"''-E"''-e"^i?^) (33) 
aoE"^ = (g-2)I]"'^-2iV<"^iV'^>^ + A^^''7V<"'^> -(5^<"(d^-r^)yl'^> 
+ e'^'^<" [ [d^ ^r^ + il^- 2A^) N^h + 2R^ ^^\] + 3H"'3 

+ {6^^°'d-,-r<" + U^°' + A<")U^^ (34) 
doN"'^ = {q5^''5 + 2^'-"s + 2e''s'^'' R^)N^^^ 

- (d^ -r^ + U^) {e'-'^^" Y.P\ ~ (5^(" R^'^ + S"^ i?^) (35) 
don = {2q-l)Q-3P~{da-2ra + 2Ua~2Aa)Q" -{T.apTl"''^) (36) 

ao«" = ^^i^ [-7i'"(ao-2g-2)r! + (G_<5";3 + 2(7-l)i;»(do-2g-2)Q^] (37) 

do^A = 2 (g + 1) r!A , (38) 



^In explicit component form these equations are available online at the URL given in Ref. 1131 . An earlier scale-invariant 
equation system (based on an orthonormal frame formulation), which employs the once-contracted second Bianchi identities 
and Weyl curvature variables, was derived by two of the authors (H.v.E. and C.U.) and given in Ref. |12|. but no specific 
choice of temporal gauge or spatial frame was introduced then. 
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where^ 

- [dp - 2rp + Up^ 2,Ap) n"'3 + e"/37 n/ . (39) 



Constraints: 



= (Cg) := 1 - f7fc - - 17 - (40) 

= (Cc)" a^iS"'' + (2,5"/3 - S"^) r'^ - S"'' - e"'''^ 7V;35 S^'' + 3Q" (41) 

= (Cj)" := (a0-r^)(iV"^ + e"'37^^)_2A;5iV"'' (42) 

= (Ca)o := (a„-2r„)r!A , (43) 



where 



(2a„ - 2r„ - 3^„) A" + \ {N^pN-^^) - ^ {N^^f . (44) 

We have also included an evolution equation and a constraint for f2A, which are a direct consequence of 
Eqs. (121, (EH and 

The role of the spatial Hubble gradient requires comment. One can use the Codacci constraint ()41|l 
to express in terms of Hubble-normalized variables. The resulting formula for involves the inverse 
of the matrix (25"^ — S"/?), and in order to avoid this algebraic complication, we propose to treat as 
a dependent variable. Choosing f — H in the commutator equations (|31|l and (|32() . and making use of 
Eqs. (|28|l and H29() . leads to both an evolution equation and a constraint for 

aor„ = [qSj + e^^^ W)rp + id^-r^ + U^){q + l) (45) 
= (C,)" := [e"^^(a3-A^)-iV"^]r^ . (46) 

These equations constitute integrability conditions for Eqs. (|28f) and H29|) . 

When we write the evolution equations and constraints as PDE by expressing do and da in terms 
of partial derivatives using Eq. H26|l . the frame components Ea^ enter into the equations as dependent 
variables. Successively choosing f = x^, i = 1,2,3, in the commutator equations (|31|l and (|32|l leads to 
an evolution equation and a constraint for -Eq,*: 

doEa' = {q6j -^J + ea^^ W)E0' -N-^daN' (47) 
= [C,,^)\fi -.^ 2[d^a-nc,~A^a)Ep^'~e^psN'^E^' . (48) 

Finally, we give the Hubble-normalized form of the gauge constraints (|20|l and (|21() : 

= {CwT ■■= [e''^''id0-rp-A0)-N^'']U^ (49) 
= {Cij)a ■■= da InAA + (r, - f/,) . (50) 

3 Gauge fixing and the Hubble-normalized state space 

In the previous section we presented a constrained system of coupled PDE that govern the evolution of 
Go cosmologies. The dependent variables are 

• the spatial frame vector field components -Eq^ 

• the spatial Hubble gradient r^, and 

• the gravitational field and matter variables S^^, A", Nap, ri, w", and Q/^. 

The system of PDE is underdetermined due to the presence of the gauge source functions 

M, N\ [/", i?" , 

which reflects the fact that there is freedom in the choice of the local coordinates and of the orthonormal 
frame. We now use this gauge freedom to specify the gauge source functions, and then proceed to describe 
some aspects of the Hubble-normalized state space. 

®We give the Hubble-normalized relativistic Euler equations, Eqs. 1361 and 1371 . in explicit form in the appendix; see 
Eqs. (T^ and IT^ . 
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3.1 Fixing the gauge 

We begin by using the coordinate freedom to set the shift vector field in Eqs. (0 and to zero: 

A^' = . (51) 

We then choose the timehke reference congruence Cq so that 

daM = . (52) 

We are then free to speciaUze the time coordinate t so that 

JV = 1 . (53) 

The effect of these choices is that the lapse function in Eq. (0 is given by = H~^, as follows from 
Eq. l|77|) . The gauge constraint l|5n|l . taken in conjunction with the above conditions, reduces to 

= {Cu)7 = (re - Ua) ^ Uc.--ra , (54) 

thus determining the frame gauge source functions Ua- The advantage of making the choices H51(l and 
(|53|l is that the temporal frame derivative do, given by Eqs. H26|l . simplifies to a partial derivative, 

do = dt. (55) 

The combined gauge choices ()51|l and H53|) have a simple geometrical interpretation in terms of the volume 
density V associated with the family of spacelike 3-surfaces S:{t — constant}, which is defined by 

V-^ := det(ea') . (56) 
Using Eq. (|51|l . the commutator equations yield 

N-'^^3, Ea'^ = - 2Aa - d,Ea' + r„ • (57) 

It follows with Eq. |H3I) that 

V^lle^^rh , (58) 

where tti = m(x') is a freely specifiable positive real-valued function of a;*, which we consider given, and 
£o is the unit of the physical dimension [length]. We thus refer to this gauge choice as the separable 
volume gauge. Note that the reduced gauge constraint H54() propagates along Eq according to Eq. (|143|l 
in the appendix. This ensures the local existence (in time) of the separable volume gauge. 
Equations (|57|l and H58(l subsequently yield the constraint 

= {Ca)^ = Aa + ^ {d,Ea' ~ ra + Ea' d^lnm) . (59) 

Finally we use a time- and space-dependent rotation of the spatial frame to relate the frame gauge source 
functions i?" to the off-diagonal components of the shear rate tensor according to^° 

(i?i, i?2, Rsf = (S23, S31, Si2)^ . (60) 

At this stage there is no freedom remaining in the choice of frame. The coordinate freedom is 

t' = t + constant , = f{x') . 

An important question, which we do not pursue at present, except in a footnote in Sec. |S1 is to what 
extent the analysis in this paper (in particular the construction of the past attractor) depends on the 
choice of temporal gauge. Here we use the separable volume gauge, as defined by Eqs. H51() and H53|) . 
which appears to be particularly well-adapted to Hubble-normalized variables. For G2 cosmologies, which 
we shall refer to later, the usual and most convenient temporal gauge is the so-called separable area gauge 
(see, e.g., Ref. J3). 

^"^In contrast to the present frame choice, one can use the frame freedom to reduce the number of variables, e.g., by 
diagonalizing the shear rate tensor. However, the present choice leads to great simplification of the equations when it comes 
to analyzing the past attractor. There are other useful choices; in particular, when one has a preferred spatial direction 
induced by an isometry. In such a case it is often advantageous to choose the iJ"-component associated with the preferred 
direction to have the opposite sign compared with the present choice. 

'^^With the exception of the special cases when the shear rate tensor is locally rotationally symmetrical or zero; when the 
frame is uniquely determined, all the Hubble-normalized variables employed are scalar invariants. 
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3.2 Hubble-normalized state space 

3.2.1 Overview 

The Hubble-normalized state vector for Go cosmologies is given by 

X ^ {E^\r^,^^p,N^p,A'',n,v'',nAf . (61) 

The evolution equations and constraints in the previous section can be written concisely in the form 

dtX = F{X,d,X,d,djX) (62) 
= C(X,a,X), (63) 

with the spatial derivatives appearing linearly (apart from the evolution equation for Tq,). A surprising 
feature of the evolution equations is that they contain second-order spatial derivatives; in this respect they 
are reminiscent of a system of quasi-linear diffusion equations. The only second-order spatial derivatives 
in the evolution equations, however, are those of the spatial Hubble gradient, didjVa, and they appear 
only in the evolution equation for itself. They arise due to the fact that in the separable volume gauge 
the deceleration parameter q contains the first spatial derivatives of r^. In fact, in the separable volume 
gauge Eq. H3()(l for q assumes the form 

Q = 2I]2 + G;^ [ (37 - 2) + (2 - 7)^2 ] 17 - - i (a„ - 2Aa) r" . (64) 

The term daq in the evolution equation for r^, Eq. H45|l . thus contains didjVa. 

A second noteworthy feature of the system of PDE H62|) is that the evolution equation for is 
homogeneous, which implies that the equation E^^ = defines an invariant set. We shall discuss the 
significance of this set later in this section. In order to clearly exhibit these aspects of the evolution 
equations, we now decompose the Hubble-normalized state vector ()61|l as follows 

X = {Ea\rc.f®Y , (65) 

where 

r = (I]„0,iV„/3,A",r!,i;",r!A)^ . (66) 
We can now write the system H62(l in a more explicit form as follows: 

dtEo.' = {qSj -J:J + ec.^''' W)Ep\ (67) 

with q given by Eq. (|64|l . and 

dtr^ - [Gj\Y) + lid^+r^)A'']rf, + Gj-'iY)df3r^-^d^{d0r^)+d^GiY) (68) 
dtVA = FA{Y)+FA^"{Y)do.YB+FA''^(Y)dc,rp + FA''{Y)ra, , (69) 

with 

da — Ea di . 

The coefhcients Gj, Gj'^, G, Fa, Fa^", Fa°''^ and Fa" are functions of the components oiY. 

3.2.2 Spatially homogeneous cosmologies 

We now discuss how the SH cosmologies are described within the Go framework. These are obtained by 
requiring that the spatial frame derivatives of the gravitational field and matter variables Y , and of the 
normalization factor H, be zero, i.e., 

daY = , r„ = . (70) 

It then follows that all the dimensional commutation functions and matter variables are constant on the 
spacelike 3-surfaces S:{t = constant}, which are thus the orbits of a three-parameter group of isometrics. 
The evolution equations and imply that the SH restrictions (fTU)) define an invariant set of the 
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full evolution equations, which we shall call the SH invariant set. Indeed, Eq. H68|l is trivially satisfied, 
and Eq. (|69|) reduces to a system of ordinary differential equations, namely 

dtYA = Fa{Y) . (71) 

The non-trivial constraints defined by (Cg), (Cc)" and (Cj)" [cf. Eqs. (^UJ^ESl] become purely alge- 
braical restrictions on Y , which we write symbolically as 

C{Y) ^ . (72) 

An important aspect of this process of specialization is that the evolution equation l|67|l for decouples 
from the evolution equation for y, which means that the dynamics of SH cosmologies can be analyzed 
using only Eqs. j7i| ) and \7S\l (cf. WE). In this context, one can think of the variables Y as defining a 
reduced Hubble-normalized state space, of finite dimension, for the SH cosmologies. 

In the SH context the restriction = defines an invariant subset, giving the so-called non-tilted 
SH cosmologies, and the Bianchi classification of the isometry group leads to a hierarchy of invariant 
subsets, some of which have been analyzed in detail in the literature. For example, the conditions 

= , A" , TV^^ = (a 7^ /3) , i?" = , = , (73) 

give the non-tilted SH perfect fluid cosmologies of class A in the canonical frame (see WE, Chap. 6, but 
with some differences in notation). 

Specializing further, by requiring the shear rate to be zero, 



, (74) 



in addition to conditions (|70|l . we obtain the FL invariant set, which describes the familiar Friedmann- 
Lemaitre cosmologies. Equations (|7U|I and (|74|l imply that = and Sap — 0, where Sap is the tracefree 
part of the 3-Ricci curvature (see App. IA.3|| . and hence that the spacelike 3-surfaces S:{t = constant} are 
of constant curvature. In addition, the electric and magnetic parts of the Weyl curvature (see App. IA.3|I 
are zero, — Sap — Ti-ap- The deceleration parameter simplifies to 

g = i (37 - 2) f7 - r^A . 

3.2.3 Silent boundary 

We noted earlier that, because the evolution equation (|67|l for is homogeneous, the equation 

Ea' = (75) 

defines an invariant set of the full evolution equations. In the introduction we discussed the notion of 
a silent initial singularity, which was introduced heuristically as an initial singularity with the property 
that the evolution along neighboring timelines decouples as the singularity is approached. In Sec. 0] we 
shall make a formal definition of a silent initial singularity, but for now we note that a key requirement 
for an initial singularity to be silent is 

lim Ea' = , (76) 

t — * — OD 

i.e., the orbit that describes the evolution of the model is past asymptotic to the invariant set Ea^ = 0. 
We will thus refer to this invariant set as the silent boundary. 

On the silent boundary, the evolution equation l|68|) for r^ simplifies to the homogeneous form 

dtra = [Gj{Y) + jraA^']rp. (77) 

It follows that the equation 

r„ = (78) 

defines an invariant subset of the silent boundary. On this invariant subset the remaining evolution 
equation H()9|) reduces to 

dtYA = Fa{Y) , (79) 

^^Noto that this does not necessarily imply limt_*_oo Ea'^ = 0. 
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which coincides with the evolution equation (|71|l for the SH cosmologies. The remaining constraints are 
purely algebraical, and can be written symbolically as 

C{Y) = . (80) 

One thus obtains a representation of the SH dynamics on the invariant set 

Ec.'=0, r„ = , (81) 

i.e., within the silent boundary. Since -Eq* = 0, however, the spatial dependence of the Hubble- normalized 
variables Y is completely unrestricted, and hence these solutions of the evolution equations and con- 
straints do not in general correspond to exact solutions of the EFE. 



4 Silent singularities and the generalized Mixmaster attractor 

In this section, we formalize the notion of a silent initial singularity, which was introduced heuristically 
in Sec.Q]^'^ We then construct an invariant set in the silent boundary that we conjecture is the local past 
attractor for Gq cosmologies with a silent initial singularity. The detailed structure of the past attractor 
in turn relies heavily on our knowledge of the asymptotic dynamics near the initial singularity in SH 
cosmologies. 



4.1 Silent initial singularities 

In terms of Hubble- normalized variables and the separable volume gauge, the spacelike initial singularity 
in a Go cosmology is approached as i — > — cxd. We now define a silent initial singularity to be one which 
satisfies 

lim Ea' = , (82) 

t — ' — OO 

lim r„ = , (83) 

t — > — oo 

and 

hm daY = , (84) 

t — > — oc 

where the E'q* are the Hubble- normalized components of the spatial frame vectors [see Eq. (jSHJ], is 
the spatial Hubble gradient [ see Eq. H29|l ] and Y represents the Hubble- normalized gravitational field 
and matter variables [see Eq. Hfififl ]. More precisely, we require that Eqs. (|82|l -184 |l are satisfied along 
typical timelines of Eq. 

One might initially think that the condition H84() is a consequence of Eq. H82() , since 

daY — Ea -Tn ■ 

However, the analysis of Gowdy solutions with so-called spikes (see Refs. ^7] and |32) shows that 
the partial derivatives dY /dx^ can diverge as t — > — oo. Thus, the requirement H84|l demands that the 
Ea' tend to zero sufficiently fast. 

We now present some evidence to justify proposing the above definition. Firstly, for SH cosmologies, 
which we have seen satisfy the restrictions (I7()|l . the evolution equation for the E^ decouples from the 
equation for Y . This evolution equation, in conjunction with the known results about the asymptotic 
behavior of the variables Y (see WE, Chaps. 5 and 6, and Ringstrom |3^), provides strong evidence that 
typical solutions satisfy^'' the remaining requirement H82() for a silent initial singularity. An example of 
an exceptional class of SH solutions, i.e., solutions for which the initial singularity is not silent, are those 
that are past asymptotic to the flat Kasner solution (the Taub form of Minkowski spacetime) , given by 

ds^ ^-dT^ + T^ dx^ + ll (d?/2 + dz^) , 

^^The concepts we propose for classifying an initial singularity as "silent" can be applied analogously to final singularities. 
^*We are indebted to Hans Ringstrom for helpful discussions on this matter. 
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where T is clock time. In terms of the dimensionless separable volume time t = lii(T/io), this line element 
reads 

ds^ = e^* (- dt^ + dx^) + dy^ + dz^ . 

from which it follows that 

lim = diag(3,0,0) . 

t — * — oo 

Secondly, further evidence is provided by recent research on G2 cosmologies, although the situation 
is clouded by the fact that an area time gauge rather than the separable volume gauge is used (but see 
the footnote in the next section about the gauge issue). Indeed, one can use the asymptotic analysis of 
vacuum orthogonally transitive G2 cosmologies (the so-called Gowdy solutions 120] ; in the present 
context see in particular Ref. |2H1) to show that the conditions (|82|) and (|84|) are satisfied along typical 
timelines, even when spikes occur. However, in general G2 cosmologies the situation is more complicated 
and further studies are needed to establish if condition (|84|l holds or if it is violated along exceptional 
timelines due to the presence of spikes. 

These results suggest that the notion of a silent initial singularity may be of importance as regards the 
description of generic spacelike initial singularities. Further support is provided by heuristic arguments of 
a physical nature, as follows. We anticipate that generic spacelike initial singularities are associated with 
increasingly strong gravitational fields, gradually approaching local curvature radii of Planck-scale order, 
which will lead to the formation of particle horizons (see, e.g., Rindler (HBjV The existence of particle 
horizons is governed by null geodesies, which satisfy 

(see Eq. (|163|) in the appendix), where are the components of the Hubble- normalized 1-forms 
associated with the orthonormal frame: 

£;",V = '5"/3- (86) 

If particle horizons form, we expect that the past-directed null geodesies emanating from a chosen point 
P will satisfy a:'(t) —> (constants) and dx^ /dt — > 0, as t — > — cxi. It follows from Eq. (|85|1 that 

dr* 
dt 

with Saph'^h^ — 1, and, hence, that 

( hm Ea"] 6" = . 

\t-*-oo J 

Since this must hold for all null geodesies emanating from P, 6" is arbitrary, implying that the limit (|82|l 
holds. In other words, we expect that the increasingly strong gravitational fields associated with a 
typical spacelike initial singularity will lead to the first condition in the proposed definition of a silent 
initial singularity. 

We now show heuristically that condition (|84|l restricts the scale of spatial inhomogeneities as the 
initial singularity is approached. If E^ tends to zero at an exponentially bounded rate as i — > —00 (as in 
SH and G2 cosmologies), the coordinate distance to the particle horizon in a direction 6" will also tend 
to zero at an exponentially bounded rate: 

Ax^ « Ea as t-^ -00 . 

The change Al^ in the Hubblc-normalized variables Y corresponding to a change Axl^ is approximated 

by 

Ai-« |^Ax^«6"a„r . 

It thus follows from the limit (|84|l that Al^ — > as t — > — cxd. In other words, the physical significance of 
the limit {84{ is that spatial inhomogeneities have super-horizon scale asymptotically as t —00, and, 
hence, up to the particle horizon scale a solution is asymptotically SH. 

^^Woei Chet Lim, private communication. 
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With the preceding discussion as motivation we now make our first conjecture. 

Conjecture 1: For almost aU cosmological solutions of Einstein's field equations, a spacelike initial 
singularity is silent. 

Proving this conjecture entails establishing the limits (|82|I - H84() . 
4.2 Stable set into the past 

We think of the evolution of the Hubble-normalized state vector X{t,x^), for fixed a;*, as being described 
by an orbit in a finite- dimensional Hubble-normalized state space. As i — > — oo, this orbit will be 
asymptotic to a past attractor, which, in accordance with the definition of a silent initial singularity [see 
Eqs. (|82|I - H84(I ], will be contained in the subset of the silent boundary defined by 

E^'^0, r„ = . (87) 

The evolution of a spatially inhomogeneous model is described by infinitely many such orbits, each of 
which is asymptotic to the past attractor. The details of the approach to the past attractor, however, 
will depend on spatial position a;% thereby reflecting the spatial inhomogeneity of the model. On the 
other hand, the evolution of an SH model will be described by a single orbit. The essential point is that 
the dynamics in the invariant set H87|l . which govern the asymptotic dynamics of both classes of models, 
is determined by the SH evolution equations and constraints, as shown in Subsec. 13.2.31 

The next step in constructing the putative past attractor is to partition the Hubble-normalized state 
vector X into stable and unstable variables, as regards evolution into the past. Firstly, within our 
framework, the BKL conjecture means that 

hm = , hm r^A = , (88) 

t — > — OG t — >~OG 

(i.e., the initial singularity is vacuum-dominated). Secondly, asymptotic analysis and numerical experi- 
ments for SH cosmologies and G2 cosmologies suggest that 

lim A" = , hm Nan (a 7^ /3) , (89) 

t — '—00 t — 00 

along a typical orbit. It is thus convenient to decompose the Hubble-normalized state vector X as follows: 

X = Xs®X^, 

where 

X,= {Ea,\ra,A'',Na.p{a^l3),n,nAf , (90) 

and 

Xu = (S„,i?",Wa,«")^ . (91) 

Here, for brevity, we have written^^ 

Xjq . — ^aa 7 Nq^ . — Nq^q^ . 

In terms of this notation, our conjectures (|82|l . 183|l . H88(l and (|89|l can be written 

lim X, = . (92) 

t — > — co 

We shall refer to the variables Xs as the stable variables, and the remaining variables X^ in Eq. 1)91(1 as 
the unstable variables. We shall provide evidence that the variables X^^ remain bounded ast—^ —00, but 
that their limits do not exist. We note in passing that further justification for the terminology "stable" 
and "unstable" in this context will be provided shortly, when we show that the variables in Xg are stable 
on the Kasner circles, while the variables in X^^ are unstable. 



^®Not to be confused with the notation used in Ref. 1251 . where Si was defined to be equal to S23, and cycle on (1, 2, 3). 
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We now list the evolution equations on the subset Xs ~ 0. Firstly, the variables Sq., R" and Na 
satisfy 

= 2(1 - E2)Ei +2(i?| - - SS"! (93) 
dtRi = [-2(1-I]2) + E2-S3]i?i (94) 
dtNi = 2(s2 + Si)iVi, (95) 

where 

S, := I iV? - I A^i (A^2 + N3)-l {N2 - N^f , (96) 
and cycle on (1, 2, 3). These variables are restricted by the Gaufi constraint H40f) . which now reads 

1 = S2 + i (Ni + Nl + Ni) -±(Ni+N2 + Ns)^ , (97) 

with 

= i (S? + + E§ + 2i?^ + 2Rl + 2Rl) . (98) 
Secondly, the evolution equation for now reads 

dtv'^ = ^ [ (37 - 4) (1 - v^) + (2 - ^){Y.p^vPv'^) ]v"-[Y."p- e"^^ (i?^ + iV^^i;^) ] , (99) 

where it is convenient to retain the index notation. We note for future use that Eq. (|99|l implies 

dtv^ = ^{l^v^)[{i^-A)v^-{^^pv''v^)] . (100) 

Although the variables are unstable into the past, it turns out that certain combinations of 
these unstable variables are in fact stable. Firstly, the limit (|92|) . in conjunction with the equation for 
dtNap {a ^ f3) and the Codacci constraint, leads to the following limits: 

lim RaNp = , P . (101) 

t — > — oo 

As a result, the subset of the Hubble-normalized state space defined by 

Xs = (102) 

is an invariant set only if the following restrictions hold: 

R^N/s = , (3 . (103) 

The essential point is that the products RaNp [a ^ (3) are stable into the past. 

Secondly, we can make use of known results about SH models to motivate another limit, in addition 
to Eq. I|1U1|I . We introduce the function 

Aat := {NiN2f + {N2Nzf + (TVsA^i)' • (104) 

If Atv 7^ 0, i.e., if more than one Na is non-zero, then Eq. 1)1031) implies i?" = 0. Then the evolution 
equations l|93|) - (|95|l reduce to the evolution equations for vacuum SH models of class A. It has been 
shown^^ that solutions of these evolution equations satisfy 

lim Aat = . (105) 

t — ^ — oo 

It is thus plausible that if Eqs. 192|) and (|101|l hold, then so docs Eq. (|105|l . We shall refer to the invariant 
set defined by 

X, = , RaNp = (a ^ /3) , A^v = , (106) 

as the stable subset into the past and make the following conjecture. 

Conjecture 2: The local past attractor A'^ for Go cosmologies with a silent initial singularity is a 
subset of the stable subset. 



Proving this conjecture entails proving the limits H92|l . (|101|l and H105|) . assuming the validity of (|82|l -184 |) . 

We believe that this conjecture can be strengthened, however. In order to do this, we need to describe 
how the Kasner vacuum solutions are represented within the present framework. 

'^^See Ringstrom 1341 for the case where the Na have the same sign (Bianchi Type— IX case). Numerical simulations 
suggest that this result is also true in the Bianchi Type- VIII case. 



4 SILENT SINGULARITIES AND THE GENERALIZED MIXM ASTER ATTRACTOR 



16 



4.3 Kasner circles 

The line element for the Kasner vacuum solutions is 

ds^ ^-dT^+ T^Pi dx^ + T^P'- dy^ + T^P"" dz^ , 
where the Kasner exponents pi , p2 and are constants that satisfy 

Pi + P2 + P3 = 1 , pl+pl+pl = i , 

and £oT is clock time. The Kasner exponents can take values that are described by the inequalities 
— 5 < Pi < < p2 < I < P3 < 1 (or permutations thereof); see Ref. (201 j P- 196. Relative to the natural 
orthonormal frame associated with this line element, the Hubble-normalized connection variables are all 
zero except for the shear rate tensor, which is diagonal and given by 

= diag (3pi - 1, 3p2 - 1, 3p3 - 1) ■ 

One can also represent the Kasner solutions relative to a spatial frame that is not Fermi-propagated, as is 
the spatial frame specified by Eq. H6U|) . Some of these alternative representations are important in what 
follows. 

Within our formulation, all possible representations of the Kasner solutions are given by 

= A" = N^p = n = v°' ^ , (107) 

r„ - , = , 6[a,Ep]' = , (108) 

with the R" given by Eq. The Gaufi constraint (gOl), together with Eqs. and imphes that 

I]2 = 1 , q = 2 . (109) 
The evolution of the non-zero variables i?Q* and Sq,^ is governed by 

dtE^' = (25„^-E«^-f e„/i?^) V (110) 
dt'E'"^ = 2eT"^<"i?-.yE'3>5 , (111) 

as follows from Eqs. l(B7|) and l(^ . 

In the physical region of the Hubble-normalized state space, i.e., det{Ea'') ^ 0, Eqs. H1U8|I imply that 
— 'Eai}{t), and that the spatial coordinate freedom can be used to obtain — Ea^{t), confirming 
that the Kasner solutions are SH and of Bianchi Type-1. On the silent boundary [Ea^ = 0), however, 
Eqs. H1U8|I become trivial, with the result that the spatial dependence of Eq,^ is unrestricted. One thus 
obtains a representation of the Kasner dynamics locally on the silent boundary, even though the line 
element, as given by Eq. H163|l in the appendix, is singular. Indeed, the Kasner dynamics on the silent 
boundary is described by the orbits that satisfy Xs — and the additional restriction Na = 0, as follows 
from Eqs. (jHOjl and l|l()7|l . We shall refer to this subset, defined by 

X, = , = , (112) 

as the Kasner set on the silent boundary. 

The evolution equations on the Kasner set are obtained by setting = 1 and S'q = in Eqs. 
and H94|l . which yields: 

dt^i = 2{Rl-Rl) (113) 
dtRi = (S2-S3)i?i, (114) 

and cycle on (1, 2, 3). Note that the evolution equations for S^t and R" decouple from that of v", discussed 
below. 

It is important to note that if the spatial frame is not Fermi-propagated (i?" ^ 0), the evolve 
in time, with E^ = 1, both on and off the silent boundary. On the other hand, if the spatial frame is 
Fermi-propagated (i?" — 0), then Eq/j is constant in time by Eq. (|lll|l . and diagonal: 



Ea/3 = diag(Ei,E2,E3) , 



(115) 
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with —2 < Si < — 1 < S2 < 1 < S3 < 2 (or permutations thereof). Thus, if the spatial frame is 
Fermi-propagated, the Kasner orbits on the silent boundary are equilibrium points of the shear evolution 
equations. Since Eq.^ is tracefree and satisfies = 1, we obtain 

Si + I]2 + E3=0, (116) 

and 

+ + . (117) 

The dynamics on the Kasner set also includes the evolution equations (|99|l for (with Naf3 — 0)- It 
follows that if the E^^ satisfy Eqs. H115() - (|117() . then the evolution equations and H100() for admit 
the equilibrium sets 

= or = 1 , 

where the latter condition is also to be supplemented with one of the six choices for v — {vi,V2,V3)'^ , 
namely 

v = ±Ea, a = 1,2, 3, (118) 

where Ei = (1,0, 0)"^, etc. Thus, there exist seven sets of equilibrium points forming circles in E^^-space, 
which we shall call Kasner circles [ the intersection of the plane H116|) with the sphere (|117|l ] , depending 
on the value of v in Eq. H118(l . which we denote by 

JC , IC±a ■ 

In addition, it follows from Eq. H99|l that for specific values of E^ subject to Eqs. 1)116(1 and H117|l . there 
are six additional lines of equilibrium points, given by 

El =37-4, vi>0 or wi < , V2^V3^0, (119) 

and cycle on (1, 2, 3), which join the various Kasner circles. 

At this stage, we digress to describe the symmetry properties of the Kasner circles. Each circle is 
divided into six equivalent sectors, which we will label according to the ordering of the diagonal shear 
components Eq, which satisfy Eqs. ((116|l and ((117|l . For example, in sector (123) these parameters satisfy 
El < E2 < E3, etc. The sectors meet at points where two of the Eq are equal. These points are of two 
types, conventionally labeled Ta (the "Taub points") and Qa, given by^^ 

Ti: (El, E2, E3) = (2,-1,-1) , Qi: (Ei, E2, E3) = (-2, 1, 1) , (120) 

and cycle on (1, 2, 3). Figure ^ represents the plane Ei + E2 + E3 = in E^-space, containing a Kasner 
circle, and showing the six sectors and the points Tq, and Qa ■ The figure also shows three additional points 
labeled Pa, which lie outside a Kasner circle, forming an equilateral triangle whose sides are tangential 
to the circle. These points, which are given by 

Pi: (Ei,E2,E3) = (-4,2,2), 

and cycle on (1, 2, 3), will be used to describe the so-called curvature transition sets. 



4.4 Transition sets 

The dynamics in the stable subset, defined by Eq. ((106|l . is essentially determined by the fact that each 
Kasner equilibrium point is a saddle point, with at least two of the nine variables Na, i?" and being 
unstable into the past. Which of these variables are unstable at a particular Kasner point can be quickly 
determined by linearizing Eqs. H95|) . H94|) and H99|) in the neighborhood of such a point. On the Kasner 
circle JC this yields 

dtNi = 2(l + Ei)iVi (121) 
dtRi = (E2-E3)i?i (122) 
dtvi = (37-4-Ei)t;i , (123) 

^*The Ta correspond to the Taub form for Minkowski spacetime in the exact Kasner solution, and the Qa correspond to 
the locally rotationally symmetrical non-flat Kasner solution. 
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Figure 1: A Kasner circle showing the six equivalent sectors and the variables that are unstable into the 
past in each sector. 

and similarly for indices 2 and 3. On the Kasner circles JC±i (where, nearby, vi = ±1 ^ Svi, Svi > 0), 
the linearized equations for A^i and Ri remain unchanged, while Eq. H99|l yields 

d^Sv, = -2^^1-^^6v, (124) 
(2-7) 

dtV2 = (Si-E2)l'2 (125) 

dtV3 = (El -Eg) 1^3, (126) 

and similarly for indices 2 and 3 on /C±2 and /C±3, respectively. It follows that 

Ni is unstable into the past 1 + Ei < 

Ri is unstable into the past <^ E2 — E3 < 

vi is unstable into the past on K. <^ 37 — 4 — Ei<0 

vi is unstable into the past on IC±i <^ 37 — 4 — Ei > 

V2 is unstable into the past on JC±i <^ Ei — E2 < 

U3 is unstable into the past on IC±i <^ Ei — E3 < , 

and similarly for indices 2 and 3. The arcs of the Kasner circles on which the variables Na and i?" are 
unstable are shown in Fig.Q] and those on which the variables are unstable are shown in Figs. [21 and |3| 
A given Kasner equilibrium point P thus has an unstable manifold (into the past).^^ The key property 
of this unstable manifold is that each of its orbits join P to some other Kasner point. The simplest such 
orbits are those on which only one of the nine variables Na, R" and is non-zero, or those on which two 
are non-zero but with extreme tilt (w^ = 1). These special orbits, which we shall refer to collectively 
as transition sets, are listed by name and symbol in Tab. 1. In this table the subscripted letter on the T 
indicates the "excited" variable. We now discuss these transition sets in turn. 

Firstly, there are the curvature transition sets T/v„ , a — 1, 2, 3, which are shown in Fig.^ In this and 
all subsequent figures, orbits are directed toward the past. For a = 1 the curvature transition orbits are 

'^^As Figs. mOl show. the unstable manifold is at most six-dimensional. 
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Figure 2: The Kasner circle /C showing the arcs on which the variables v"' are unstable into the past. 
Figure 2a shows the case | < 7 < 2 and figure 2b the case 1 < 7 < |. 



(213) 




Figure 3: The Kasner circles IC±i showing the arcs on which the variables are unstable into the past. 
The variable vi is unstable on the boldface arc to the right of the line Sn =87 — 4. 

given by 

(l-fc)(2-E2) = (l + fc)(2-S3) , (127) 

where /c is a parameter that satisfies — 1 < fc < 1. This relation follows from Eqs. (|93|l and H95|l with 
= A^2 = ^3 = R"- In the spatially homogeneous setting these orbits describe the Taub vacuum Bianchi 
Type-II solutions, and determine the past attractor for vacuum and non-tilted SH models of Bianchi 
Type-VIII and Type-IX (in a Fermi-propagated frame; see WE, Fig. 6.6 on p. 138, and p. 143-7, but 
note differences in labeling) . These curvature transitions link different "Kasner epochs" , according to a 
transition law for the Kasner exponents first found by BKL (see Ref. p. 535-7; also WE, p. 236). We 
derive this transition law in App. I A. 41 

Secondly, there are the frame transition sets, Tfi^ , which are shown in Fig. [3 For a — 1 they are 
given by 

Si = fc , 

where fc is a parameter that satisfies — 2 < fc < 2. In the spatially homogeneous setting these transition 
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Table 1: The transition sets 



Name Symbol Transitions 



curvature 




/C -> /C (w" = 0) 
or 


1) 


frame 




/C ^ /C (w" = 0) 
or 

IC±a IC±a (f ^ = 


1) 


tilt 




K- IC±a or vice 
(Sa/j fixed) 


versa 


extreme-tilt 




IC±a IC±i3 

(Sa/3 fixed, = 1) 



























Figure 4: The curvature transition set TJvi- The sets T/Va and T/Vg are obtained by cycling on (1, 2, 3). 

sets map a Kasner solution into a physically equivalent Kasner solution through rotation of the spatial 
frame by tt/2 about one of its axes. For example, the transition sets Tjjj result in the interchange 
S2 < S3 — > E3 < as can be seen by comparing Figs. ^ and 




Figure 5: The frame transitions sets Tr^ (5a), Tr^ (5b) and Tr^ (5c). 

Thirdly, there are the tilt transition sets, 7^^, which are shown in Fig. El They are simply lines with 
Sq, constant and one of the v" non-zero. Whether the orbits join a point on /C to a point on one of the 
Kasner circles IC±a, or vice versa, depends on the values of the (see Figs. |21 and The reversal of 
direction of these orbits is governed by the six lines of equilibrium points given by Eq. (|119|l . 
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Figure 6: The tilt transition set % 



Finally, there are the extreme-tilt transition sets. Let us consider the subset — 1 (with fixed Sq). 
Then = e", where e" is a unit vector, which may be parametrized according to 

e" = (cos sin 1? cos 93, sin I? sin <^)^ , < I? < TT , < <y9 < 27r . (128) 

We then obtain from Eq. (|99|l a simple dynamical system for the polar angles (p}, given by 



= — a^sin2i? := ^ {— pi + P2 cos^ ip + siii^ ip) 

dt(p = sin 2^3 P := -§(^2-^3) 



It is easily seen that if Si < S2 < S3 (or equivalently pi < P2 < Ps), then the past attractor of this 
dynamical system is given by ip} — {^tt, ^n} and {i?, ip} = {^tt, |7r}, and, hence, 

lim e" = ±£;3 ; (130) 

t — ^—00 

similarly for other orderings of S^. Thus, extreme-tilt transition sets are orbits that lie on the extreme-tilt 
sphere = 1 in t;"-space, with = = R", and Eq, fixed. Figure[3shows the directions corresponding 
to El < E2 < S3, i.e., the arc (123) on the Kasner circles. The other cases can be obtained by 
interchanging 1, 2 and 3. The points B±a in Fig. [7| correspond to the points on the Kasner circles IC±a 
determined by the values of the S^. The directions of the orbits joining the points B±a depend on the 
ordering of the Eq,. 



4.5 Structure of the past attractor 

We are now in a position to state our conjecture concerning the local past attractor of the Go evolution 
equations for generic initial data satisfying the constraints. We introduce the following notation: 

Tn—HTn^ , T]i=UTfi^ , Ttilt=LI%^ , "T^xtromo = U 7^^„^ , /Coxtrcmc = U /C±q . (131) 

Thus, T)v is the union of all curvature transition sets, Tr is the union of all frame transition sets, etc. 
(see Tab. 1 for the complete list of transition sets). 

We now make the following conjecture concerning the local past attractor . 

Conjecture 3: The local past attractor for Go cosmologies with a silent initial singularity is 
given by 

A — ICU /Ccxtromc U T/v U Tr U Ttut U Extreme ■ (132) 
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Figure 7: The extreme-tilt transition sets on the extreme-tilt sphere = 1, = = R", with fixed 
and Si < E2 < S3. 

The essential property of the various transition sets is that they define so-called infinite heteroclinic 
sequences on the past attractor, i.e., infinite sequences of Kasner equilibrium points joined by transition 
sets, oriented into the past. In particular, a typical Kasner point will be the starting point for infinitely 
many heteroclinic sequences (infinitely many, because at least two transition sets (three, when 7 < |) 
emanate from each Kasner point). The significance of the heteroclinic sequences is that an orbit that 
is asymptotic to the past attractor (for given values of the a;*) will shadow a heteroclinic sequence, and 
hence the cosmological model will be approximated locally by a sequence of Kasner states. 

The conjectured past attractor, A~ , is a proper subset of the stable subset, defined by Eq. H106|l . 
Thus, in addition to the stable variables H9()|l. various expressions involving the unstable variables Na, 
R" and will tend to zero on the attractor, even though the limits of these variables as t ^ —00 do 
not exist. The desired expressions depend on the following quantities: 

and 

An := iNiN2f + (iVzTVa)^ + (NsN^)^ , 
together with analogous quantities and A„ for i?" and w". The past attractor is then characterized 



by the following hmits [in addition to Eqs. (|H21)-(|H2l and 

limt_,_oo N^R^ = (133) 

limt^_oo (1 - v^) (N^ -f i?2) = (134) 

limt^_oo (Ajv,Afl) = (0,0) (135) 

limt^^o, (A„(iV2+i?2), A,(l-i;2)) = (0,0) (136) 

limt^_oo V1V2V3 = . (137) 



The limits ()133ll and (|134|l imply that if (i?" , respectively) is active, than R" {Na , respectively) must 
be close to zero, and likewise either t;^ or 1 — t;^ must be close to zero, as on the T/v„ {Tu^ ) transition sets, 
respectively. The limit (|134|l also implies that if is not close to or 1, then Na and i?" must be close 
to zero, as on the tilt transition sets Ty^ . The limit (|135|l implies that at most one Nq. and at most one 
R" can be active simultaneously, as on the TJv^ and Tr^ transition sets. The limit l|136(l implies that if 
two are active simultaneously, (i.e., A^, 7^ 0), then Na and i?" must be close to zero and v"^ must be 
close to 1, as on the extreme-tilt transition sets Ty^y^. 

The conjectured structure of the past attractor in terms of Kasner equilibrium points and transi- 
tion sets, as given by Eq. (|132|l . or as described by the limits (|Hlll-l|Hl|), (ESl and (|133l) - (|137|l . embodies 
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the notion that as one approaches the attractor into the past along an orbit, the probabihty that more 
than one of the nine unstable variables Na , i?" and v" is active during any one transition (except for the 
extreme-tilt transition sets on which two are non-zero, but are constrained by = 1) tends to zero. 
Evidence that the probability of multiple transitions involving pairs such as (-Ri, iVi) or (i?i, i?2) tends to 
zero as t — > — oo is provided by numerical simulations for non-tilted SH models of Bianchi Type-VII-^yg, 
reported by Hewitt et al [^ . 

We shall refer to the local past attractor , defined by Eq. H132|) . as the generalized Mixmaster 
attractor, since it generalizes the past attractor for the so-called Mixmaster models (SH models of Bianchi 
Type-IX; see WE, p. 146, and Tab. 2 to follow), making precise the heuristic notion of an oscillatory 
singularity. 

5 Cosmologies with isometries 

In Sec. ^ we proposed a detailed structure for the past attractor for Gq cosmologies with a silent initial 
singularity [ see Eq. (|132|l ] . Classes of cosmologies that admit an isometry of some sort are described by 
invariant sets of the Hubble-normalized state space. In this section we exploit this fact to predict the 
structure of the past attractor for these more specialized models, thereby providing a link to much recent 
research. 

For models with an isometry it is possible that one or more of the nine unstable variables Na, i?" or 
v" is required to be zero, leading to two possibilities. 

(i) The initial singularity is oscillatory . 

This possibility occurs if each arc of the various Kasner circles has at least one unstable variable 
(refer to Figs.mSJ. The attractor will then include all available transition sets, and the evolution 
into the past along a typical timeline will be described by an infinite sequence of Kasner states, 
possibly of a more specialized nature than for Gq cosmologies. 

(ii) The initial singularity is Kasner-like. 

This possibility occurs if at least one arc on one of the Kasner circles has no unstable variables. 
The arc(s) in question then form the past attractor, and the evolution into the past along a typical 
timeline will be described by a specific Kasner state. A cosniological solution with this type of 
singularity is also referred to as being asymptotically velocity (-term) -dominated, a term that has its 
origins in the work of Eardley, Liang and Sachs 10 and Isenberg and Moncrief 26 . 

We now present various classes of cosmologies with isometries, whose initial singularities have been 
discussed in the literature, and give the conjectured past attractor in terms of our formulation. The 
specific nature of the isometry determines whether the initial singularity is oscillatory or Kasner-like. 

Firstly, we consider SH cosmologies, which, as we have seen in Subsec. 15.2.21 can be described by the 
finite-dimensional reduced state space defined by the variables Y . Since the definition of the generalized 
Mixmaster attractor A~ involves only the variables Y , the attractor also exists as an invariant set in 
this reduced state space. Indeed, we conjecture that in this context A~ is the past attractor for the 
general class of SH cosmologies, and that it will thus contain the past attractors for the three special 
classes of SH cosmologies with an oscillatory singularity that have been analyzed in detail to date. In 
Tab. 2 we give these three special classes of SH cosmologies and list the key references. These papers 
use Hubble-normalized variables, but there are some differences in the labeling of variables compared to 
the present paper. It is notoriously difficult to prove rigorous results about oscillatory singularities and 
little progress has been made until recently, when Ringstrom |34j , in a remarkable piece of mathematical 
analysis, rigorously established the existence of the past attractor for the class of non-tilted Bianchi 
Type-IX cosmologies, by proving the required limits, conjectured earlier by WE (see p. 146-7). 

Secondly we consider spatially inhomogeneous cosmologies. Most recent research on the initial singu- 
larity has been restricted in two ways: 

(i) the spacetime is assumed to have compact spatial sections, 

(ii) the energy-momentum-stress tensor is assumed to be zero (vacuum solutions). 



5 COSMOLOGIES WITH ISOMETRIES 

Table 2: Perfect fluid SH cosmologies with oscillatory initial singularity 



24 



Class of models Non-zero Past attractor 

unstable variables 



non-tilted Type-VIII and Type-lX TVi , iVj, iVg /C U Tat 

(WE, p. 146, Ringstrom |34)) 

non-tilted Type-VF ^/g TVs, i?3 /C U TAfg U Tfl^ U Tr^ 

(Hewitt et al |2S1) 

tilted Type-II A^s, i?3, fi ^ U /Ci U T^Vg U Ti?^ U U T^,, 
(Hewitt et al pi)) 



The first restriction is made because it enables one to prove results about the global existence of solutions, 
ft does not, however, affect the structure of the past attractor, since it is determined by the dynamics along 
individual timelines. In view of the BKL conjecture, namely that matter is not significant dynamically as 
the initial singularity is approached, one might believe that the second restriction can be made without 
loss of generality when determining the past attractor. This conclusion is not valid, however. Our analysis 
leads to the conjecture that the past attractor for vacuum Go models is in fact the much simpler set given 

by 

A-,, ^ICUTnUTr, (138) 

since the Hubble-normalized state vector for vacuum models does not contain the peculiar velocity variable 
v" , which implies that the extreme Kasner circles /Coxtrcmc and the transition sets Tint and T^xticmc cannot 
be part of the past attractor [see Eq. H132(l ]. Nevertheless, determining the vacuum past attractor is an 
important first step in determining the attractor for non-vacuum models. 

In Tab. 3 we list the classes of vacuum spatially inhomogeneous cosmologies whose initial singularity 
has been studied. In each case we can predict immediately whether the singularity will be oscillatory 
or Kasner-like. In the table we give the past attractor for each class, which is a subset of the general 
vacuum attractor A^^^^.'^'^ 

One class of vacuum models is not included in Tab. 3, the so-called polarized C/(l)-symmetric mod- 
els (Berger and Moncrief These are Gi cosmologies for which the single spacelike Killing vector 
field is hypersurface-orthogonal. The reason for this exclusion is that the dynamical consequences of 
the hypersurface-orthogonality condition are not compatible with our choice of spatial frame, given by 
Eq. Hfi()|l . These models could be incorporated by making a different choice of spatial frame, as discussed 
in footnote 10, but we will not pursue this matter here. 

It should be noted that in the papers listed in Tab. 3 the conclusions about the dynamics near the 
initial singularity are not expressed in terms of a past attractor: we have reformulated their results within 
our dynamical systems framework, and at this stage most of the results about the past attractor have 
not been rigorously established. The papers referred to use a metric-based approach^^ instead of the 
orthonormal frame approach. Some of them make use of the so-called Fuchsian algorithm to establish the 
asymptotics at a Kasner-like initial singularity (see Refs. [2H| and [23) while others rely on a Hamiltonian 

^"At this stage the reader might be concerned with the fact that these models have not been studied in the separable 
volume gauge. However, we believe that, due to asymptotic silence, our discussion is "gauge robust", i.e, that the local 
asymptotic temporal behavior is not affected by the choice of temporal gauge. To make this more substantial we note that 
the choice Af = 1 and = was not necessary for obtaining our picture of the past attractor. Any sufficiently smooth 
choice J\f{X) such that J\f is positive and bounded on the attractor does not change the flow on the past attractor and thus 
one would obtain the same results as the choice J\f = 1 yields; can be similarly generalized. We also note that these 
are not necessary conditions, and that even wider sets of gauge choices are allowed if one takes into account the detailed 
structure of the EFE. 

■^^In a recent paper 0, however, the Gowdy models are analyzed using scale-invariant variables introduced in Refs. 1231 
and 
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Table 3: Past attractor for vacuum spatially inhomogeneous cosmologies with isometrics. In the first three 
cases, the initial singularity is Kasner-like and in the last two cases, the initial singularity is oscillatory. 



Class of models Non-zero Past vacuum 

unstable variables attractor 

Polarized Gowdy all zero K. 
= diagonal G2 



(Isenberg and Moncrief [211 ) 

Unpolarized Gowdy 
= OT G2 
(Kichenassamy and Rendall |28| 

Polarized T2-symmetric 
= G2 with one HO KVF 
(Isenberg and Kichenassamy j27| 'l 

T2-symmetric 
= generic G2 
(Berger et al ^) 

U ( 1 )-symmetric 
= generic Gi 
(Berger and Moncrief [7]) 



Unpolarized Gowdy A^3,-Ri arc(T2(3i) C /C 

= OT G2 



Polarized T2-symmetric R3 arc(T3Q2T'iQ3) C /C 

= G2 with one HO KVF 



T2-symmetric A^s, i?3 JCU T/v., U Tr^ U Tr^ 

generic G2 



[/(l)-symmetric all non-zero /C U T/v U Tr 

= generic Gi 



formalism and the so-called method of consistent potentials to predict whether the initial singularity will 
be oscillatory or not (see Refs. |H] and 0). In this approach, the transitions between Kasner states 
are described heuristically as bounces off potential walls determined by the Hamiltonian. Some of these 
papers also describe numerical simulations that display a finite number of Mixmaster oscillations. 



6 Concluding remarks 

In this paper we have developed a mathematical framework for analyzing the dynamics of Go cosmologies, 
and in particular the BKL conjecture discussed in the Introduction. A key step was the introduction 
of scale-invariant variables, using the Hubble scalar defined by a timelike reference congruence as the 
normalization factor. One of the principal advantages of Hubble-normalization lies in the behavior of the 
dynamical variables as the initial singularity is approached: the dimensional variables diverge, while, for 
at least a generic family of solutions, the Hubble- normalized variables remain bounded. 

The structure of the evolution equations and constraints led to the introduction of the silent boundary 
in the Hubble-normalized state space and enabled us to define a silent initial singularity. The next step 
was to construct the generalized Mixmaster attractor^ which makes precise the notion of an oscillatory 
initial singularity in a Go cosmology, while having a simple geometrical structure (see Figs. and 
allowed three precise conjectures on early cosmological dynamics to be formulated (Conjecture 1 in 
Subsec. im Conjecture 2 in Subsec. W?R and Conjecture 3 in Subsec. r4.5|l . The construction of the past 
attractor also highlights and clarifies the important role of SH dynamics in the Gq context. Indeed, there 
is now considerable evidence, both numerical and analytical, that SH dynamics influences the asymptotic 
dynamics of spatially inhomogeneous cosmologies near the initial singularity in a significant way. Our 
formulation places this relationship on a sound footing: the local past attractor for Go cosmologies with 
a silent initial singularity is the past attractor for SH cosmologies . We are now in a position to restate 
the BKL conjecture in a precise form: 
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For almost all cosmological solutions of Einstein's field equations, a spacelike initial singularity 
is silent, vacuum- dominated and oscillatory. 

Proving this conjecture entails establishing all the limits associated with Conjectures 1, 2, and 3 in Sec. 01 
As a first step, one would have to complete the proof for the SH models, begun by Ringstrom [311 • A 
natural second step would be to consider the simplest class of spatially inhomogeneous models with an 
oscillatory initial singularity, namely the generic G2 models (see Tab. 3), restricting consideration to 
vacuum solutions for simplicity. Analyzing the role of the spatial derivatives in a neighborhood of the 
silent boundary will be a major step in this analysis, and will clearly present a formidable challenge. 

This unifying statement incorporates certain fundamental physical ideas about singularities, partially 
supported by known examples and theorems. It is useful to revisit the conjectured physical behavior in 
a way that highlights various aspects of the situation: 

(i) The generic cosmological initial singularity is a strong-gravity phenomenon, and so should be linked 
to trapped surfaces, which intuitively capture the notion of ultra-strong gravitational fields (and 
thus also to the standard singularity theorems). 

(ii) The generic cosmological initial singularity is likely to be a spacelike curvature singularity because a 
null singularity will be very special and timelike singularities will by their nature intersect relatively 
few worldlines of matter (but confirming this will depend on implementing a good measure on the 
space of cosmological models, which is needed in any case in order to put on a firm footing all talk 
about probabilities). 

(iii) The generic spacelike curvature singularity is a scalar curvature singularity, since non-scalar curva- 
ture singularities require fine tuning of initial data. 

(iv) If the energy conditions are strictly obeyed, the curvature singularity is generically Weyl curvature 
dominated, at least when vorticity in the matter fluid is not significant (this is not the case if the 
energy conditions are just marginally satisfied, as exemplified by stiff perfect fiuids, but these are 
not physically likely states). 

(v) The strong-gravity regime associated with the initial state leads to particle horizons, and spatial 
inhomogeneities are constrained to have super-horizon scale as the initial singularity is approached. 

(vi) Increased strength of the gravitational field and the collapse of particle horizons lead to asymptotic 
silence, and on the scale of the particle horizon solutions therefore are asymptotically SH. 

(vii) The past attractor describing asymptotic spatially inhomogeneous dynamics is thus given by the 
generalized Mixmaster attractor. 

(viii) Final singularities are in essence the time reverse of initial singularities, and so we expect the above 
ideas to apply there too, and conversely that ideas from gravitational collapse can throw some 
light on cosmological initial singularities. In particular, in generic gravitational collapse, angular 
momentum plays a fundamental physical role, so this should also be true in many time-reversed 
cases, i.e., at the initial singularity; but when this is true, matter is dynamically important, in 
contrast to the cases considered above. 

Each of these issues needs to be investigated and given a precise mathematical statement; e.g, the last 
may relate to the hypothesis that the tilt transition sets can be interpreted as physical/dynamical effects 
of vorticity in the matter fluid and associated transverse peculiar velocity components. In each case 
we wish to link our results and conjectures to physical ideas. The challenge is to explain the difference 
between the past attractor for vacuum Go models and the past attractor for perfect fluid Go models 
[compare Eqs. H132I) and (|138|l ] using physical principles. It is often stated that "matter (energy) does 
not matter" in the approach to the initial singularity — this view is reflected in our past asymptotic 
limits for the matter variables Q and Qa. But, perhaps — and this is rather heuristic and speculative at 
this stage — , "matter linear momentum does matter" and/or "matter angular momentum does matter". 

In the end the major physical statements are 
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• Ultra-strong gravitational fields will occur in the early Universe, associated with local restrictions 



• Propagating gravitational iua,ves are not im,portant in the cosmological context, hut tidal forces are, 
and indeed are often more important than the gravitational fields caused directly by the matter. 

• The relation between tidal forces and vorticity in the matter fluid is unclear and may contain some 
of the most interesting physics. 

The relation between them is that — if our conjectures are correct — in the early Universe, energy and 
information mainly propagate along timelike worldlines rather than along null rays. When matter moves 
relative to the irrotational timclikc reference congruence, as must be the case when vorticity in the matter 
fluid is dynamically important, then the energy and information will flow with the matter. The primary 
effect of the gravitational field is in determining the motion of the matter through Coulomb-like effects; 
on the other hand, the effect of the matter on the gravitational field is primarily through concentrating 
that field into small regions, while conserving the constraints which embody the Gaufi law underlying the 
Coulomb-like behavior. The effect of spatial curvature is to generate oscillatory behavior in tidal forces 
as this concentration takes place, as seems to be characteristic of generic cosmological initial singularities; 
but this is not wavelike in the sense of conveying information to different regions, it is just a localized 
oscillation. 

It is issues such as these that need to be investigated when further developing the themes studied 
here. 

Acknowledgments 

We thank Woei Chet Lim for many helpful discussions and for performing numerical simulations, and 
Joshua Horwood for preparing the figures. C.U. and J.W. appreciate numerous stimulating discussions 
during February 2003 with participants of the workshop on "Mathematical Aspects of General Relativity" , 
which helped clarify our understanding of Go cosmologies and contributed to the final form of this 
paper. This workshop was held at the Mathematisches Forschungsinstitut, Oberwolfach, Germany, whose 
generous hospitality is gratefully acknowledged. H.v.E. acknowledges repeated kind hospitality by the 
Department of Physics, University of Karlstad, Sweden. C.U. was in part supported by the Swedish 
Science Council. J.W. was in part supported by a grant from the Natural Sciences and Engineering 
Research Council of Canada. 

A Appendix 

A.l Propagation of constraints 

Propagation of dimensional constraints: 



on causality. 



eo[(C'com)a/3(/)] 



(C'com)a/3[eo(/) ] - 2H {C co-m) afiif) + '^'^[0^ (Ccom) l3]f{f) 



eoiCc) 
eoiCcr 



- 2e^5[„f^T(Ccom)/3]*(/) (139) 

- 2H (Cg) + 2 (e„ + 2u^ - 2a„) {Cc)" + e"^"" (Ccom)a/3(fi^) (140) 

- [ 4HS''0 + a'^p - e"^0 fi^ ] {Ccf - ^ {S"^ - 2w") (Cg) 



+ \ e^^^ [ep +up- 3a^) (Cj)^ - | n"^ {C,f + i n/ (Cj)" (141) 



eo(Cj) 



- [ 2i? <5"/3 - a'^p - W ] {C,f - \ e"t'-< (C™,n),3-, {H) 



+ \ e"''^ (Ccon.)/35(a^') - \ e^^' {C.om)M{^"s) " (Qom)"/3(0^) • (142) 



Propagation of dimensionless gauge fixing condition: 



(143) 
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A. 2 Hubble-normalized relativistic Euler equations 

Upon substitution of the matter variables (|25(l . Eqs. H36|l and l|37|l assume the exphcit form 

dM = --^t;"a„f} + G;M2G+<7-(37-2)-(2-7)t.2-7(E„^t,"«'3) 

- 7 (a^ - 2r„ + 2Ua -2Aa)v" +j In G+ ] f7 (144) 

dov" = - yf^ dpv" + (5"^ d;; In G+ - ^ (1 ~ v^) 5°"^ {dp In 17 - 2rp) 

7 

+ G-J [ (7 - 1) (1 - v^) {dpyP) - (2 - 7) Bp In G+ 

+ (2 - 7) (1 - w^) (5/5 In r! - 2r/3) + (87 - 4) (1 - i;^) 

7 

+ (2 - 7) {^^yv'') + G_ (;7/3i>'3) + [G+ - 2(7 - 1)] {Apv^) ] 

- I]"/3 + e"/37 i?'^ - f/" - A" + e"'^'' Nps . (145) 

Using Va = ve", BqC" = 1, we easily obtain from Eq. 11451) 

(jT — 

- Ilull (1-^2) (a„ In f7 - 2r„) 

7 

+ (37 - 4) «2 - (S„;3f %'^) - G_ [Uc^v^] - 2 (7 - 1) (A^t,") ] (146) 

doe" = - z; d^e" + - In G+ - - ^ (1 - f') p"^ (da In » - 2rff ) 

w u 7 

- p'^fi + s"/3 R'^ -\ P^P -V p-^n Af^ +v s^'p N>^^ . (147) 
Here p"-p :— S"p — e"ep and :— e"^^ e'^. 

A. 3 Hubble-normalized curvature variables 

The Hubble-normalized 3-Ricci curvature of a spacelike 3-surface S:{t = constant} is defined through the 
symmetric-tracefree and trace parts. The trace part ftk was given in Eq. 1)44(1 . while the tracefree part is 
given by 

Sa.p - i e^'(a (^hl - " ^^hl) Np)S + i (^(a " r(„) A^) + | 7V(„^ iV^)^ - i TV/ iV(„/5) . (148) 

The quantities Sap and satisfy the Hubble-normalized twice-contracted 3-Bianchi identity given by 

0={dp- 2rp - 3Ap) S"" - e^/'T Nps 5/ + i S""" (dp - 2rp) ■ (149) 
Employing Eq. p48|) . we can write the evolution equation for "Eap in the alternative form 

aoS"'' = {q^2)^"f^ ^3{S"^ -n"'^) + e''^<°'[2R-yY.'^\~N^\Us] 

+ ((5''<"a-y-r<"-f C/<" + A<")C/^> . (150) 

The conformal curvature properties of a spacelike 3-surface S:{t — constant} are encoded in the Hubble- 
normalized 3-Cotton-York tensor 

Ca^p := e^\^ (9,^1 - 2r|^| - A\^^)Sp^s - SN^^^Sp^^ + N^^ S^p . (151) 

The Hubble-normalized Weyl curvature variables take the explicit form 

£af3 — Sap + ^'Sa0 ~ ^^{a'^ ^p)'f — ^^aP (152) 

Hap = le'^\a{d\j\ - ri^i - A\^i)j:p)s - N(^a'^ Y.p)^ + ^ Nj'^ T,ap , (153) 
with Sap defined in Eq. H148|l . 
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A. 4 Curvature transitions 

Although the relation H127|l implicitly gives the rule for the relationship between two Kasner epochs, 
k is not particularly suitable for explicitly describing the "Kasner transformation law" for curvature 
transitions. However, that law can be elegantly obtained in the present dynamical systems framework as 
follows. The solutions on the TJvi -subset are determined by 



2 



and 

Si = -4 + 3Z, S2 = 2-|r+Z, Eg = 2 - § r_ Z , 9tZ = - 2(1 - E^) Z 



wherer± := IzbVl - a^, a e [0,1], is a constant, and where l-E^ | [(q!-2)Z + 2] [(q: + 2)Z-2]. An 
orbit starts (with time direction reversed toward the past) from a Kasner point where Z = Z_ =2/ (2+q:), 
and ends at a Kasner point where Z = Z^ = 2/(2 — a). That is, is a parameter on the individual 
orbits that increases monotonically from Z^ to Z+ toward the past, while a labels the different orbits. 
The value Z = determines the point P\ outside the Kasner circle in Fig. It is possible to express 
the constant a, and, hence, r-t, in terms of the standard Kasner parameter u > 1 (see, e.g., BKL 
p. 528), where we assume that we are considering orbits that originate from sector (123) and thus that 
Pi < Vi < P3- Then (r_|_,Q;,r_) = 2(u^,m, 1)/(1 + u^), and thus a might be viewed as a compactified 
Kasner parameter. The above formulae directly yield the standard transformation laws from the initial 
Kasner exponents {pa} to the final Kasner exponents {p^ } for orbits originating from sector (123): 

Curvature transitions from other sectors than (123) are easily obtained through permutations. 
A. 5 Inverting the commutator equations 

The commutator equations can be used to solve for the connection variables in terms of the frame 
variables. Let us introduce the dual frame variables e°'i of Cq.*, and their Hubble-normalized counterparts 
E°'i := H e^i. The dual frame variables can be used to conveniently describe the line element, giving 

ds^ = H-^ [ - JV^ dt^ + E", Ef^j {NWt + da;^ {N^dt + d.T^ ) ] . (155) 

As Ea" and E"'^ satisfy the relations E"'^ Efj' = S^p and E", E^^ = (5,,^ one finds with Eq. gT)! that the 
E°'i evolve according to 

daE''^ = - {qS'^p - S"^ - e''^f3W)E^, +M-^E''jd^N' . (156) 

On the other hand, from appropriately inverting Eqs. (|47|) . (|48ll and H5U|) we obtain the explicit expressions 

(157) 
(158) 
(159) 
(160) 
(161) 
(162) 

When working in the separable volume gauge, determined by conditions (|53|l and (|51|l . the line element 
takes the form 

ds^ =H-^[-Ae + 5apE'',Ef^jAx'dx^] . (163) 

Combining the above result for A^ + with the constraint (|59|) makes it then possible to solve for 
and Aa separately: 



1 


= \E^,doEa' + \N-^d,N' 




= -6^l^aE\doEp)' - N-H^^aE\dp)N' 


i?" 




Aa + Ta 


= i E^i daEp"^ — i d-iEa^ 








= da InTV . 



To 



i E^i EJ djEp' + ^Ea'^^\nm (164) 



3 ^ » '^J^P I 3 



Aa = I E^, Ea' djEp' - i d^Ea' - i Ea' ^^\nm. (165) 
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